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I. INTRODUCTION 

In recent papers the general formalism to study poly¬ 
tropes for anisotropic matter has been presented, both 
in the Newtonian [T] and the general relativistic regimes 
[113]. The motivations to undertake such a task were, 
on the one hand, the fact that the polytropic equations 
of state allow to deal with a variety of fundamental as¬ 
trophysical problems (see Refs. I3H3I] and references 
therein). 

On the other hand, the local anisotropy of pressure 
may be caused by a large variety of physical phenomena 
of the kind we expect to find in compact objects (see Refs. 
[22H32j , and references therein for an extensive discussion 
on this point). 

Among all possible sources of anisotropy (see [22] for a 
comprehensive discussion on this point), let us mention 
two which might be particularly related to our primary 
interest. The first one is the intense magnetic field ob¬ 
served in compact objects such as white dwarfs, neutron 
stars, or magnetized strange quark stars (see, for exam¬ 
ple, Refs. [53H37] and references therein). Indeed, it is 
a well-established fact that a magnetic field acting on a 
Fermi gas produces pressure anisotropy (see Refs. [38l - 
SI] and references therein). In some way, the magnetic 
field can be addressed as a fluid anisotropy (however, it 
should be observed that for the thermodynamic pressure 
of degenerate particles under the conditions of Landau 
quantization, it remains isotropic |42] 1. 

Another source of anisotropy expected to be present 
in neutron stars and, in general, in highly dense matter, 
is the viscosity (see [iSHSO] and references therein). At 
this point it is worth noticing that we are not concerned 
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by how small the resulting anisotropy produced by the 
viscosity might be, since as we shall see below, the occur¬ 
rence of crackings may happen even for slight deviations 
from isotropy. 

An alternative approach to anisotropy comes from ki¬ 
netic theory using the spherically symmetric Einstein- 
Vlasov equations, which admits a very rich class of static 
solutions, none of them isotropic (Ref. [5TH53] and ref¬ 
erences therein). The advantages or disadvantages of ei¬ 
ther approach being related to the specific problem under 
consideration. 

The theory of polytropes is based on the polytropic 
equation of state, which in the Newtonian case reads 

P = Kpl = Kp\^^l\ ( 1 ) 

where P and po denote the isotropic pressure and the 
mass (baryonic) density, respectively. Constants K, 7 , 
and n are usually called the polytropic constant, poly¬ 
tropic exponent, and polytropic index, respectively. 

In the general relativistic, anisotropic, case, two possi¬ 
ble extensions of the above equation of state are possible, 
namely: 

1 . 



(2) 


(3) 


where Pr and p denote the radial pressure (see below), 
and the energy density respectively. 

As it should be expected, the assumption of either (H 
or ([^ is not enough to integrate completely the field 
equations, since the appearance of two principal stresses 
(instead of one as in the isotropic case) leads to a system 
of two equations for three unknown functions. 
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Thus in order to integrate the obtained system of equa¬ 
tions we need to provide further information about the 
anisotropy, inherent to the problem under consideration. 
For doing that, in [5] a particular ansatz has been as¬ 
sumed, which allows for specihc modelling. This method 
links the obtained models continually with the isotropic 
case, thereby allowing us to bring out the influence of 
anisotropy on the structure of the object, even for very 
small anistropies. Here we shall adopt the above men¬ 
tioned ansatz. However it should be stressed that our 
models are presented with the sole purpose of illustrat¬ 
ing the occurrence of crackings; the natural way to ob¬ 
tain models, consists of providing the specific informa¬ 
tion about the kind of anisotropy present in each specihc 
problem. 

The concept of cracking was introduced to describe the 
behavior of a huid distribution just after its departure 
from equilibrium, when total nonvanishing radial forces 
of different signs appear within the system [55] . 

Thus, if we take a snapshot of the system, just af¬ 
ter leaving equilibrium (on a time scale smaller than the 
hydrostatic time scale), we say that there is a cracking 
whenever this radial force is directed inwards in the inner 
part of the sphere and reverses its sign beyond some value 
of the radial coordinate. In the opposite case, when the 
force is directed outwards in the inner part and changes 
sign in the outer part, we say that there is an overturn¬ 
ing. Further developments on this issue may be found 
in [55H55] . As should be clear at this point, the concept 
of cracking is closely related to the problem of structure 
formation [HIHHH] . 

In the example examined in [55] it appears that crack¬ 
ing occurs only in the locally anisotropic case whereas in 
the perfect fluid case, once out of equilibrium, the conhg- 
uration tends either to collapse or to expand as a whole. 
Furthermore it has been shown in EH, that for a wide 
range of models, cracking appears only if, in the process 
of perturbation leading to departure from equilibrium, 
the local anisotropy is perturbed. This result suggests 
that fluctuations of local anisotropy may be the crucial 
factor in the occurrence of cracking. Thus, in the case 
of an initially locally isotropic configuration, the appear¬ 
ance of cracking shows that even small deviations from 
local isotropy may lead to drastic changes in the evolu¬ 
tion of the system as compared with the purely locally 
isotropic case m- 

It is the purpose of this work to study in some detail the 
conditions under which, cracking appears in anisotropic 
polytropes. We shall consider the two possible cases of 
anisotropic polytropes mentioned above, and our models 
are the same as those considered in [5]. 

At this point, we have to make some comments about 
the perturbative scheme, leading to departures from equi¬ 
librium. For the cracking to occur, it is imperative that 
the perturbations introduced into the system take it out 
of equilibrium. In other words we shall consider, exclu¬ 
sively, perturbations under which the system is dynami¬ 
cally unstable. One way to assure this, is to assume that 


the value of the ratio of specific heats of the fluid is not 
equal to the critical value required for marginal (neutral) 
dynamical stability. So, under perturbations, the config¬ 
uration either collapses or expands. Another completely 
equivalent way to achieve this instability consists in as¬ 
suming that, under perturbations of density and local 
anisotropy, the radial pressure of the system maintains 
the same radial dependence it had in equilibrium. It is 
obvious that in this case the hydrostatic equation will no 
longer be satisfied. 

It should be clear that assuming this lack of response 
of the fluid, i.e. the inability to adapt its radial pressure 
to the perturbed situation, is equivalent to assuming that 
the pressure-density relation (the ratio of specific heats) 
never reaches the value required for neutral equilibrium. 

In this work, we shall consider two possible ways to per¬ 
turb the energy density and the local anisotropy, namely: 

• We shall perturb the parameter K and the param¬ 
eter describing the anisotropy. 

or 

• We shall perturb the parameter n and the param¬ 
eter describing the anisotropy. 

As we shall see, cracking occurs for a wide range of values 
of the parameters. 


II. RELEVANT EQUATIONS AND 
CONVENTIONS 

A. The field equations 

Even though polytropes are static configurations, we 
have to keep in mind that the appearance of cracking 
occurs when the system is taken out of equilibrium. Ac¬ 
cordingly, we shall consider spherically symmetric distri¬ 
butions of collapsing fluid, which we assume to be locally 
anisotropic, bounded by a spherical surface E. The static 
limit is trivially recovered from the equations below. 

The line element is given in Schwarzschild-like coordi¬ 
nates by 


= e''df - e^dr^ - {dO^ + sin^ edcj)^) , (4) 

where v[t^ r) and A(t, r) are functions of their arguments. 
We number the coordinates: = t; = r; x^ = 

9; x^ = (j). 

which in our case read |55| : 


-SttFO 




(5) 
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- SttTI = -SttTI 


2 A + A(A-z>) 



2v" + 


- \'v' + 2 


/-y 


, (7) 


-S^Toi = - ( 8 ) 

r 

where dots and primes stand for partial differentiation 
with respect to t and r respectively. 

In order to give physical significance to the Tj^ compo¬ 
nents we apply the Bondi approach m- 
Thus, following Bondi, let us introduce purely locally 
Minkowski coordinates (r, x, y, z) 


Note that the coordinate velocity in the (t, r, 0, (j>) system, 
dr/dt, is related to w by 




dt 


Feeding back (10 131 into 
tions in the form 


(14) 

we get the field equa- 


p + PrUj'^ __1 f ^ 

I — Uj2 J.2 ® I ^2 j. 


(15) 


Pr + ^ f ^ 

1 — oj^ ^ r 


(16) 


dr = e'^^^dt, dx = e^^^dr, dy = rd9, dz = r sin 9dcj). 

Then, denoting the Minkowski components of the energy 
tensor by a bar, we have 

^0 ^0 'T^l 'T^l rTi2 rTi2 

^0 — ^0 ^ -^1—^1^ ^2 — ^2 ^ 

T| = T|, Toi = e-(-^+^'>/^Toi. (9) 




2A-f A(A-!>) 


o-A r 


2i/" + y" - AV + 2^^^^ 


(17) 


(p -f _ A 

1 — r ' 


(18) 


Next, we suppose that when viewed by an observer mov¬ 
ing relative to these coordinates with proper velocity oj in 
the radial direction, the physical content of space consists 
of an anisotropic fluid of energy density p, radial pressure 
Pr and tangential pressure P±. Thus, when viewed by 
this moving observer the covariant tensor in Minkowski 
coordinates is 


/ p 0 0 0 \ 

0 P,. 0 0 1 

0 0 Pj_ 0 ■ 

\ 0 0 0 P± / 

Then a Lorentz transformation readily shows that 


rpO _ ^0 _ 

-^0 — -^0 — 


= 


Tqi = 


1-W2 ’ 

(10) 

Pr + pU)^ 

1-W2 ’ 

(11) 

1 

II 

Fs" 

II 

(12) 

(p -1- Pr)oje^''^^P'^ 



toi — 


At the outside of the fluid distribution, the spacetime 
is that of Schwarzschild, given by 


ds^ = ^1 — dt^ — , ^^2M\ ~ ■*" 9d(l)^) . 

(19) 

In order to match smoothly the two metrics above on the 
boundary surface r = we require the continuity of 

the first and the second fundamental forms, across that 
surface. 

These last conditions imply 


= 1 - 


= 1 - 


[Prh = 0 , 


2M . . 

1 

rs 

(20) 

2M 

•> 

ty. 

(21) 

0 , 

(22) 


where, from now on, subscript E indicates that the quan¬ 
tity is evaluated at the boundary surface E. Eqs. (201, 
(21) and (22) are the necessary and sufficient conditions 
for a smooth matching of the two metrics (|^ and (19) 
on E. 

The energy momentum tensor may be written as: 


Tfj. 1 , = {p+ P±) UfrUr, - Pj_gtJ.u + {Pr “ P±) (23) 
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where denotes the four velocity of the fluid and is 
a unit spacelike vector, radially directed. These vectors 
are defined by 


which follows from (15) and (16) in the static case, we 
may write 


, _(. "■ +4 . ^ ^ ^ (29) 

r (r — Ira) r 


,-v/1 


„-A /2 


(l-a;2)^/2’ (1-^2) 


2 \ 1/2 


, 0,0 


(24) 


where the mass function in(r), as usually, is defined by 

e~^ = 1 — 2mlr, m(r) = An f pr^dr. (30) 

Jo 


=- 1-/2 


= -A/2 


(l-w2)l/2’ (1_ ^2)1/2 


0 , 0 


(25) 


Next, it will be useful to calculate the radial component 
of the “conservation law” 


=0 


(26) 


Polytropes are static fluid configurations, which satisfy 
either of the equation ([^ or (|^ . The full set of equations 
describing the structure of these self-gravitating objects, 
in both cases, were derived and discussed in [5]. 

All the models have to satisfy physical requirements 
such as: 


P > 0, — < 1, — < 1. (31) 

P P 

In what follows, we shall very briefly review the main 
equations corresponding to each case. 


which in the static case becomes 


i? = p; + ^(p + p,)-^^^^—^=0, (27) 

representing the generalization of the Tolman- 
Oppenheimer-Volkof equation for anisotropic fluids 
[22] . Alternatively, using 

_ m + AnPrT^ 

2 r(r-2m) ’ ^ ^ 


B. Case I 

Assuming Eq. (§, let us introduce the following vari¬ 
ables 

a = Prc/pc, r = ^IA, = Anpc/a{n + l), (32) 

^0 = Po/Poc, = w(r)A3/(47rpc), (33) 

where subscript c indicates that the quantity is evaluated 
at the center. At the boundary surface r = (^ = ^s) 

we have = 0. 

Then, the generalized Tolman-Opphenheimer-Volkoff 
equation becomes 




1 — 2(n -|- l)av/^ 


(1 — na) + (n + 1)q!4'o 




2A4'o”^ r l-2a{n + l)v/^ ' 

Prc{n -1-1) (1 — no) + {n + l)Q!'ko 


(34) 


where A = P± — P^. 

On the other hand we obtain from the mass function 


In this case, conditions (31) read 


definition (30) and Eq. (15) 


m! = Anr^ p, 


(35) 


or 


dv 


^ 2 ^n(l _ „Q, _|_ 


(36) 


na < 1, 


a4'o 

1 — na + naTo 


< 1 , 


3u/C^ -f Q;4'g+^ 
'kg (1 — na -I- naA>o) 


(37) 
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C. Case II 

In this case the assumed equation of state is ^ , then, 
introducing 


the generalized Tolman-Opphenheimer-Volkoff equation 
becomes 


= P/Pc. 


(38) 


, 2 ^ 


1 — 2(n + l)av/^ 


I + a^' 


+ V + - 


Prcin + 1 ) 


1 — 2a(n + l)v/^ 
1 + 


= 0 , 


(39) 


r 


and from Eq. (351 




In this case, conditions (311 read: 

3v 

p > 0, < 1, ^ 

^ ’ £3x]/n 


- 1 < 1 . 


(40) 


(41) 


Equations (34), (361 or (39), (40), form a system of two 


first order ordinary differential equations for the three un¬ 
known functions: 4'(4'o), r:. A, depending on a duplet of 
parameters n, a. Thus, it is obvious that in order to pro¬ 
ceed further with the modeling of a compact object, we 
need to provide additional information. Such informa¬ 
tion, of course, depends on the specific physical problem 
under consideration. Here, we shall further assume the 
equation of state used in [U . 


III. PERTURBING THE ANISOTROPIC 
POLYTROPE 

Let us now consider an anisotropic polytrope satisfy¬ 
ing the equation of hydrostatic equilibrium ( |27[ ). Fur¬ 
thermore, our distribution satisfies the equation of state 
proposed in [^, i.e. 


A — Ci^P + Pr) 


m -I- AnPrV'^ 
r{r — 2m) 

where C is a constant, producing 


(42) 


R=^ + h 

dr 


m -I- 47rr^Pr 


r[r — 2m) 


{p + Pr) — 0, 


(43) 


with 


h = l-2C. 


(44) 


In order to observe an eventual cracking, we have to 
take the system out of equilibrium. For doing so we shall 


perturb the energy density and the anisotropy, through 
the parameters K and h, or, alternatively, through the 
parameters n and h. Below we describe this procedure 
for the cases I and II. 

A. Perturbation in the case I 

In this case the following equations apply 

= and p = po+nKpl'^^/'^. (45) 

Then, perturbing the energy density and the local 
anisotropy via K and h, 


it follows that 


K —> K = K + 5K, 
h —> h = h + Sh, 


Pr = = ^Pr, 


p = po + 'nkp]^^l'^ = Po + nfiPr, 


with 




K 


(46) 

(47) 

(48) 

(49) 

(50) 


where the tilde denotes the perturbed quantity. 
Introducing the dimensionless variable 


Airpi 


(51) 


we obtain from (45)-(50) 


A = + f [(1 - »a) + + 1)*„| 

(52) 

From the above, it follows that up to first order, we 
may write 
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SR = R{^, l + SP, h + 5h,v + 5v) = R{^, l,h^v) + 



P = 1 ■'^^ + 
h = h 



V = V 


p = i 
h = h 

V = V 



P = 1 -Sh. (53) 
h = h 

V = V 


where 


R{i,l,h,v) = 0 , 


as it must be, since the system is at equilibrium in the 
unperturbed state. 

(54) Thus, we may write 


5R = 




V = V 


P = 1 + 

h = h 

V = V 



p = 1 

h = h 

V = V 


(55) 


Then, using (52) we obtain 


im 


^ = 1 
h = h 

V = V 


(dR\ 


(1 — na) + a(n + l)'I'o 

ydv j 

P = l 

u u 

(1 — 2Q;(n + l)u/^)2 


{n + l)(u + 2^3 c^4'”+^) + - no) 

1 — 2a;(n + l)u/^ 


[l + 2a\n + l)e'^l+P 


(56) 


(57) 


V = V 


Im 


vl/^ 


/(= 1 = ^ 1(1 -"»)++ i)»ol ■ 

h = h 


V = V 


(58) 


r 


Furthermore, we have 


producing 


1 ri 

V=— (POc + nPPrc'i'o)dl 

Pc Jo 


and 


-(I 


.sp, 


(59) 


Sv = na[ / Sp = naF{p)Sp, 


with 


F(C) = / f4/^+ide 

(60) From the equations above, we may write 


(61) 


( 62 ) 




d^ 

d^ + 


ha 

^2(1 - 2a{n + l)u/^) 


[(n + l)(u + 2e^a4'((+i)^o + " na) 
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(1 — na) + a{n + l)4'o 
1 — 2a{n + l)v/^ 


+ ^ [(1 - na) + a{n + l)^'o] 


[l + 2a2(n + l)^2^JJ+ij 

'u + ^3^«'o+^ 


6P 


I — 2a{n + 


5h. 


(63) 


r 


It would be more convenient to use the variable a;, defined in terms of which, (63) becomes 
by 


a;=^, A = r^A = ^^, 


(64) 




ha 


A { dx x{Ax — 2a{n + l)u) 
1 — na + a(n + l)4^o 


[(n + l)(u + 2A^a;^a4'Q+^)4'oA^a;^4'o+^(l - na) 


Ax 


Ax — 2a{n + l)t> 


(1 + 2a^(n + l)y4^a;^4'”+^)nFi 


6P 


^ (1 — na + a(n + l)'I'o) 




Ax — 2a(n + l)n 


6h. 


(65) 


r 


Now, for the cracking to occur, once the system has been 
perturbed, there must be a change of sign in 6R. More 
specifically, it should be positive in the inner regions, and 
negative at the outer ones (in the inverse case an over¬ 
turning is produced). This last condition implies that 
6R = 0 for some value of ^ in the interval [0 ,^e], imply¬ 
ing in turn: 


where 


r = 


Sh = —r(5/3. 


dR j_ OR Tp( 


( 66 ) 


dv ' 
dR 
dh 


/3 = 1 

V = V 

h — h 


(67) 


Another perturbative scheme may be proposed, based 
on the perturbation of the energy density and the 
anisotropy through the parameters n and h. In this case, 
following the same line of arguments as in the previous 
scheme, we have 


Pr = ( 68 ) 

p = Po + nPr = Pock'd + (69) 


n —> h = n + Sn 

h —> h = h + Sh. 


(70) 

(71) 


Then, assuming that the radial pressure remains un¬ 
changed under perturbation (what ensures the departure 
from equilibrium) 


P — P 

± r — J- Tt 


(72) 


we have 


P, = P, = 


P = POc'^"" +hPr. 


(73) 


(74) 


Then, replacing into the hydrostatic equlibrium equa¬ 
tion, it follows that 






h 

e 


[(1 - na)4'(( + {h + l)aTo+^] 


V + C^a4'o+^ 

(1 — 2a(n -I- l)u/^) 


(75) 





























Proceeding exactly as in the previous scheme, we obtain 


5R = 


+ 




(1 - 2a(n + l)u/C) 


° [1 — na + (n + 1 )q;'I'o] 


{v + - na) In('I'o) + a'I'o) 

V + 


,n+l 

-0 


1 — 2a(n + l)u/^ 


6h, 


1 — na + (1 + n)a'I'o 
1 — 2a(n + l)i'/C 


(l + 2a"f (n+l)^'”+i)i^i„ 


or, 


SR = 




Ax{Ax — 2a(n + l)n) 

1 — na + a(n + l)'I'o 


{v + yl^a;^ad>Q''‘^)((l — na) In('I'o) + a^'o) 


Ax 


Ax — 2a(n + l)n 


{1 + 2a^{n + l)A^x^’i>^+^)Fin 


6n 


+ (1 — na + a(n + l)'I'o) 

Ax 


V + A^x^ad'g^^ 

Ax — 2a(n + l)v_ 


5K 


(77) 


with 


with 


Fin= / e2d'^[(l-na)ln(«'o)+a«'o]de (78) 

Jo 

Again, for the cracking to occur we must have 6R = 0 
for some value of ^ in the interval [0 ,^e], implying : 


with 


r = 


Sh = —T5n, 


dR I dR rp 

jMyAUhhLif 

OR 

dh 


n = n 

V = V 

h = h. 


B. Perturbation in the case II 


For this kind of polytropes, we have 


Pr = Kp^ + ^/^, p = 


Po 


(1 + ATpi/")”' 

Then, we shall write for the radial pressure 


(79) 


(80) 


(81) 


G= ( ^ 

Pc 


1 + 1 /n 


(83) 


thereby ensuring that the radial dependence of Pr re¬ 
mains unchanged after perturbation. 

We shall now proceed to perturb our polytrope, first 
through the parameters K and h. Then, following the 
same line of arguments as for the case I, we have 


Pr = kpl+^^'-G = PPr, 


P = 


Po 


Po 


(i+iFp^”)" [i+pKpy^Y 


or, using 


p = p + 5p = p + 
/3 = 1 -f 5j3, 


dp) 0= 


SP 


0=1 


we obtain 


(84) 


(85) 


( 86 ) 

(87) 


F, = Kpl+^/^G, 


(82) 


p = p + nPr{P - 1) 


( 88 ) 
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Feeding back the above expressions into (431 we may Carrying out a procedure similar to the one described 
write in the previous subsection we obtain 


(89) 


vk” f d’F 
SR=^\^ 


ha 


x4 ( dec x{Ax — 2a(n + l)r;) 


[4'(u(n + 1) + A^ec^4'"(l + a{n + 2)4')) 


+ Ax ( ^ 1 X ) (1 + + l)Q!^A^a;^4'"+^)nF2 

Ax — 2a{n + l)u 


\I/" 

SP + ^(1 + Q;4') 

J\.X 


V + 


_Ax — 2a{n + l)r;)_ 


Sh. (90) 


where (64) has been used, and 


F2= (91) 

Jo 

If the perturbation is introduced via the parameter n, 
then the result is: 


SR = 




x{Ax — 2a{n + l)r;) 


{v + A'^ec'^a4'”+ ) ln(4') + Ax 


1 + 0:4' 


Ax — 2a{n + l)z 


(l + 2(n+l)aMV^”+^)Fi 


6n 


V]/" 

+ ^(l+a*) 


V + j4^a;^a4'”+^ 
Ax — 2a(n + l)r;) 


(92) 


with 


F2n = / ln(4')^"de 


(93) 


IV. RESULTS 


We shall now proceed to apply the formalism sketched 
in the previous section, to detect and analyze the occur¬ 
rence of cracking in the family of polytropes defined by 
the equation of state (42). We shall consider both types 
of polytropes (1 and 11)7 both schemes of perturba¬ 
tion (through K and h, and through n and h). 

To do so, using a four order Runge-Kutta method, we 
shall integrate (34), (36), and (39), (40), for any duplet 


of parameters n, a and n, h , for which there exists a 
boundary surface. Such duplets have been determined in 
[5] (see figures (3, 4) for the case 1, and (8, 9) for the case 
11 , in that reference). Also, the absolute numerical values 


of the perturbations of the parameters were taken of the 
order of 10“®, the results are not qualitatively sensitive, 
to changes in some orders of magnitude, below or above, 
that figure. Finally, the integrals [7^ |M| were 
numerically evaluated by means of tl 
for example [68]). 

Next, the obtained 4', 4'o (which of course coincide 
with the results in [5]) are used to evaluate ( j^ and 
( |77| ) for the case I, and (|^ and (j^ for the case II. We 


le trapezoid rule (see 


have considered the whole range of values of the param¬ 
eters leading to bounded configurations, whose physical 
variables satisfy the physical requirements (31). 


We shall now analyze the obtained results for each case. 


A. Cracking for polytropes of type I 

Figures (1-7) summarize the main results obtained for 
the polytropes of type 1. The first four figures (figs. (1- 
4)) correspond to the case where the perturbation is car¬ 
ried on through the parameter AT, whereas the remain- 
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FIG. 1. SR as a function of x for n = 1, h = 1.5 and F = 1.6. 
Curves a, b, c and d correspond to a = .83, .85, .87 and .90, 
respectively. 



FIG. 2. SR as a function of x for n = 1, h — .5 and 
F = 1. Curves a, b, c, d and e correspond to a = 
.80, .83, .85, .87 and .90, respectively. 

ing three figures of this case (figs.(5-7)) describe the be¬ 
haviour of the system which has been perturbed through 
the parameter n. 

In figure 1 we observe the cracking fo /i = 1.5 and the 
different values of a indicated in the figure legend. It ap¬ 
pears that the strongest, and deepest, crackings are asso¬ 
ciated with the largest values of a. We also observe that 
in the outer regions, the models exhibit overturnings. In 
these latter cases, smallest values of a are associated with 
strongest overturnings. 

Figure 2 depicts a similar situation, but for h = .5. 
However in this case crackings are stronger and overturn¬ 
ings are weaker. 



FIG. 3. SR as a function of x for n = 1.5, a = .390, F = .6 
and different values of h. Curves a, b, c, d, e and f correspond 
lo h — .5, .7, .9, 1.1, 1.3, and 1.5, respectively. 

The situation described in figures 1 and 2 is representa¬ 
tive for a wide range of parameteres (for which there exist 
bounded configurations satisfying the required physical 
conditions). 

Figure 3 exhibits further the close relationship, exist¬ 
ing between the occurrence of cracking and the type of 
anisotropy. Indeed, the smallest values of h correspond 
to the deepest and strongest crackings, and to the weak¬ 
est overturnings. Furthermore, for the largest values of 
h (curves e and f) there are not crackings at all, only 
overturnings. 

Finally figure 4 exhibits the occurrence of cracking for 
conhgurations that are locally isotropic before perturba¬ 
tion, thereby illustrating the fact that even slight devia¬ 
tions from local isotropy, may be sufficient to produce a 
cracking. 

Figures 5-7, which describe the results of perturbation 
through the parameter n, lead to conclusions qualita¬ 
tively similar to the ones mentioned above. The strongest 
and deepest crackings, and the weakest overturnings, are 
associated with largest values of a 


B. Cracking for polytropes of type II 

The results of this case are depicted by figures (8-14). 
The first four of these, correspond to the perturbative 
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FIG. 4. 5R as a function of x for n = 1, h = 1 and 
r = 1.2. Curves a, b, c, d and e correspond to a = 
.80, .83, .85, .87 and .90, respectively. 



FIG. 5. 5R as a function of x for n = 1, h — 1.5 and F = 
1.4. Curves a, b and c correspond to a = .85, .87 and .90, 
respectively. 


scheme through the parameter K whereas the last three 
correspond to perturbations of n. 

Figure 8 describes the situation, for different values of 
h. For the choice of parameters in this figure, we ob¬ 
serve that cracking only appears for the smallest value of 
h, whereas for the two largest values, there is only over¬ 
turning. For h = .7;.9 neither cracking nor overturning 
occurs. 

This behaviour is conhrmed in figure 9. Here we see 
how the strongest and deepest crackings are associated 
with the smallest values of h. For some minimal value of 
h, and larger, only overturnings are observed. These are 
strongest for the largest values of h. 

Figures 10 and 11 also confirm the role of h in the 
occurrence (or not) of crackings. In figure 10, for h = .5, 
there are only crackings, which are strongest and deepest 
for the largest values of a. In contrast, hgure 11 only 
depicts overturnings for h = 1.5 and different values of 
a. 

Figures 12, 13 and 14 describe the situation resulting 
from the perturbation of n. However, unlike the previ¬ 
ous cases, for these three figures we have assumed neg¬ 
ative values of F. This explains why they only exhibit 
overturnings, and why these overturnings are stronger 
for smaller values of h (just the opposite of the situation 
observed for positive values of F). 


V. CONCLUSIONS 

We have investigated the conditions under which, gen¬ 
eral relativistic polytropes for anisotropic matter, exhibit 
cracking and/or overturning, when submitted to fluctu¬ 
ations of energy density and anisotropy. 

Thus, we have shown that cracking and/or overturn¬ 
ing occur for a wide range of the parameters. For both 
types of polytropes, the main conclusions are basically 
the same, namely: the strongest and deepest crackings 
occur for the smallest values of h and greatest values of a. 
Also, the strongest overturnings result from the largest 
values of h. 

A distinct feature of all the models studied here, is the 
fact that the “core” (the inner part), and the “envelope” 
(the outer part), respond differently to different degrees 
of anisotropy (different h). This fact was already pointed 
out in HUBS]. 

It is important to stress that the occurrence of a crack¬ 
ing, has direct implications on the structure and evolu¬ 
tion of the compact object, only at time scales that are 
smaller than, or at most, equal to, the hydrostatic time 
scale. This is so because as already mentioned, what we 
do is to take a “snapshot” just after the system leaves 
the equilibrium. To find out whether or not the system 
will return to the state of equilibrium afterward, is out 
of the scope of our analysis, and would require an inte- 
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FIG. 6. 5R as a function of x for n = 1, /i = 1 and 
r = .8. Curves a, b, c, d and e correspond to a = 
.80, .83, .85, .87 and .90, respectively. 

gration of the evolution equations for a finite period of 
time, greater than the hydrostatic time. 

However, all this having been said, it is clear that the 
occurrence of cracking would drastically affect the future 
structure and evolution of the compact object. 

Accordingly, we would like to conclude this work by 
speculating about possible scenarios where the occur¬ 
rence of cracking might be invoked, in order to under¬ 
stand the related observational data. 

One of these situations could be the collapse of a su- 
permassive star. The occurrence of cracking at the inner 
core, would certainly change (in some cases probably en¬ 
hance) the conditions for the ejection of the outer man¬ 
tle in a supernova event. This will be so for both the 
“prompt ” [TDiin] and the “long term” mechanisms [Ta¬ 
ils]. 

Also, one is tempted to invoke cracking as the possible 
origin of quakes in neutron stars [7BH75| . In fact, large 
scale crust cracking in neutron stars and their relevance 


in the occurrence of glitches and bursts of x-rays and 
gamma rays have been considered in detail by Ruderman 
(see m and references therein). 

Evidently, the characteristics of these quakes, and 
those of the ensuing glitches, would strongly depend on 
the depth at which the cracking occurs. In this respect 
it is worth noticing, the already mentioned fact, that the 
depth at which the cracking may appear, is highly depen¬ 



FIG. 7. SR as a function of x for n = 1, h = .5 and 
r = .8. Gurves a, b, c, d and e correspond to a = 
.80, .83, .85, .87 and .90, respectively. 

dent on the parameter h, which measures the anisotropy 
of the pressure. 

However we would like to emphasize that our aim here 
is not to model in detail any of the scenarios mentioned 
above, but just to call attention to the possible occur¬ 
rence of cracking in such important configurations, as 
those satisfying a polytropic equation of state, and its re¬ 
lationship with fluctuations of local anisotropy. In other 
words, whatever the origin of the anisotropy would be, or 
how small it is, cracking may occur, which would drasti¬ 
cally affect the outcome of the evolution of the system. 
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